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O ' I. INVARIANT AND NON-INVARIANT GAUSSIAN RANDOM MATRIX ENSEMBLES 

(N 

Any random matrix ensemble (RME) is determined through the probability distribution function (PDF) P(M) 
that depends on the matrix entries M nm . An important special class of random matrix ensembles is given by the 
PDF which is invariant under rotation of basis M — * TMT -1 : 
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P(M) cx exp [-Tr V(M)] , (1) 



where V(M) is and arbitrary function of M analytic at M = 0. The invariance is ensured by the trace Tr in front of 
the matrix function V(M). The RME with the PDF of the form Eq.{T]) will be referred to as invariant RME. 

In general the PDF Eq.{T]) corresponds to a non-trivial correlation between fluctuating matrix entries. However, 
^ ■ there is one extremely important case when 

• ' (i) the matrix M is Hermitean M = H = and 
Si (ii) V(H)=aH 2 . 
d In this case 
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so that all matrix entries fluctuate independently around zero. This is the celebrated Gaussian random matrix 
ensemble of Wigner and Dyson (WD). 

Note that Gaussian random matrix ensembles can be also non-invariant. The generic non-invariant Gaussian RME 
■ is determined by the PDF of the form: 
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(2) 



In this case each matrix entry fluctuate independently of the other but with the variance which depends on the indices 
n, m that label the matrix entry. The simplest Gaussian non-invariant ensemble is the Rosenzweig- Porter ensemble 2 
for which 

_ f a, n^m ( ) 

Anm ~\Aan = m {6) 

It is remarkable that both the classic WD ensemble and the Rosenzweig- Porter ensemble allow for an exact solution-^. 

Physically, an invariant random matrix ensemble describes extended (but phase-randomized) states, where the 
localization effects are negligible. In contrast to that any non-invariant ensemble accounts for a sort of structure of 
eigenfunctions (e.g. localization) in a given basis which may be not the case in a different rotated basis (remember 
about the extended states in the tight-binding model which are the linear combinations of states localized at a given 
site). 

In particular the problem of localization in a quasi- 1 wire can mapped onto the non-invariant banded RME with 
the variance matrix equal to: 

Aim =cxp \-\n-m\/B] (4) 



This model can be efficiently mapped onto a nonlinear supcrsymmetric sigma model and solved by the transfer matrix 
method^. 
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Finally, we mention a critical power-law banded random matrix ensemble (CPLB-RME) for which the variance 
matrix is of the Lorcnzian form: 
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(5) 



This model (not yet solved) possesses a fascinating property of multifractality and is an extremely accurate model for 
describing the critical states at the Anderson localization transition point in dimensionality d > 2. 

II. PARAMETRIZATION IN TERMS OF EIGENVALUES AND EIGENVECTORS 



Consider a Hermitean N x Amatrix H = EF. The physical meaning is mostly contained in the eigenvalues E n 
of this matrix and also in the unitary matrix U = (U^) _1 whose n-th column is an n-th normalized eigenvector 
^„ = {^E'„(r)}. Therefore it is sensible to parametrize the matrix H in the following way: 

H = UELP\ (6) 

where E = diag{£"„}. Then instead of N(N + l)/2 independent entries of the Hermitean matrix H one will deal with 
N(N — l)/2 independent variables of the unitary matrix U plus N eigenvalues. 

For invariant ensembles the PDF is independent of the eigenvector degrees of freedom and is determined only by 
the eigenvalues. In particular for the classic WD ensemble it reduces to: 



P(H) cx exp 



(7) 



However the change of variables involves also computing the Jacobian of the transformation Eq. (|164p . The easiest 
way of computing it is to compute the form: 



Tr(dH) 2 = Tr (2df Edf E - 2E 2 (df) 2 + (rfE) 2 ) = 2 £ (E n - E m ) 2 \df nm \ 2 + £(d£„) s 



(8) 



where 



df = U 1 dV = -ciU t U 



The set of f nm , (n > m), are N(N — l)/2 natural "coordinates" related to eigenvectors. Then the Jacobian of the 
transformation dH — > di dE is given by 
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yD, H is real 
D, H is complex 
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A 2 = [] (E n - E m f 



(9) 



III. JOINT PROBABILITY DISTRIBUTION 



According to Eqs . (J7J) , (J9j) the entire joint probability distribution function of eigenvalues and eigenvectors for an 
arbitrary invariant ensemble Eq.([T]) takes the form: 



rfHP(H) = dtd{E n } exp 
where A is the Vandcrmond determinant: 
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(11) 



A remarkable property of this distribution is that it is independent of the eigenvectors and depends only on eigenvalues. 

Since the Vandcrmond determinant is vanishing if any two eigenvalues coincide E n = E m (two columns of the 
determinant are equal) the coincidence of two eigenvalues is statistically improbable. This is the basic property of 
the random matrix theory which is called level repulsion. 
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IV. LEVEL REPULSION: POOR MAN DERIVATION 

In order to understand the physical origin of level repulsion let us consider a situation where occasionally two levels 
are very close to each other \E\ — E 2 \ <C A, where A is the mean level separation. Then it is enough to consider only 
one block of the random matrix: 

si V 
V* e 2 

The true energy levels of this two-level system are well known: 

E ia = £ -±±^± 1 -^{e l -e 2 y + \V\\ (12) 

The two-level correlation function which is the probability density to find a level at a distance ui from the given 
one, is given by: 

R{u) = J de 1 de 2 VV [s(oj - y/( £l - e 2 ) 2 + \V\ 2 ) - S(u - £1 + e 2 )] P(ei, e 2 , V), (13) 

where 

VV-{ dV ' /?=1 (14) 

Small energy difference \E\ — E 2 \ <C A implies that both e% — e 2 and \V\ are small. Then the PDF P(si,e 2 , V) can 
be considered independent of e% — e 2 and \V\. Thus integrating the (5-functions over si — e 2 we arrive at: 



R(ui) = / VV 



9(u>-\V\). (15) 



Apparently this integral is convergent and the power counting immediately leads to: 

J: 2 ^ 

The simple analysis above illustrates two important points. One -physical- is that the level repulsion is nothing 
but the avoided level crossing which is well known in quantum mechanics and which is caused by the | | 2 term in the 
square root in Eq. (fT2|) . The other one - formal- is that the pseudo-gap in R{uS) near uj = is the effect of the phase 
volume T>V and the power of uj depends on the number of independent components of V which is 1 in case V is a 
real number and 2 if V = 3tV + iQV is a complex number. 

It is known that the algebra of real and complex numbers allows only one further step of generalization. This is 
the algebra of quaternions: 

TQ = 1, Tl,2,3 = ^1,2,3) (17) 

where ci,2,3 are 2x2 Pauli matrices. It appears that this generalization makes sense in the context of random matrices 
too. Namely, one can consider random matrices which entries are real quaternions, i.e.: 

3 

V = ^iTi, (18) 
with real components £j. This generalization corresponds to (3 = 4 in Eqs. (|32l39[) . 

V. TIME-REVERSAL SYMMETRY AND THE DYSON SYMMETRY CLASSES 

It turns out the the parameter (3 is related with the time-reversal symmetry. In order to see this we note that the 
time-reversal operator T should obey a basic property 



T 2 = al, \a\ = l 
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(time reversal applied twice leaves the wave function unchanged). As the time reversal operator should involve the 
complex conjugation of wave function one may write: 

T = KC 1 (19) 
where C is the complex conjugation operator and K is an operator such that 

KCKC = KK* =al. 

But K must be a unitary operator (as the norm of the wave function must be conserved). That is why 

K* K T = 1. 

From these two conditions one finds: 

K = aK T = a(aK T f = a 2 K. (20) 

Thus we conclude that 

a 2 = 1 =>a = ±l. (21) 

For spinless particles (or particles with even spin) we have: 

T 2 = 1, (22) 

and K can always be chosen to be a unity operator K = 1. 
However, for particles with half-integer spin 

T 2 = -1, (23) 
and K is not an identity operator. In particular for spin-^ particles: 

The physical meaning of this operator is very simple: it flips the spinor. 
The time-reversal symmetry 

= H 

in the cases Ed.(|22p and Eq. (f2"3")l implies, respectively 

H = H* (25) 

and 

H = -KH*K. (26) 

In the first case time reversal symmetry requires the Hamiltonian matrix H to be real, which corresponds to /3 = 1. In 
the second case one can do a simple algebra exercise and show that the condition Eq. (j2l))) is fulfilled if the Hamiltonian 
matrix H has entries of the form Eq. (|18|) with real coefficients £j . As was already mentioned this case corresponds to 
= 4. 

Thus we see that the case (3 — 1 (known as the orthogonal ensemble) corresponds to the particles with an even 
spin and a Hamiltonian that preserves the time reversal symmetry. The case (3 = 4 (known as the symplectic 
ense?7i6/e)corresponds to particles with an odd spin and a spin- dependent Hamiltonian that preserves time reversal 
symmetry. In particular it applies to a system with the spin-orbit interaction. The case (3 = 2 (known as the unitary 
ensemble)does not assume any definite relationship between H and H*, and thus the time-reversal symmetry must 
be broken (e.g. by magnetic field or magnetic impurities). 

According to the initial idea of Wigner and Dyson all systems with complex interactions do not possess any 
symmetry but possibly time-reversal symmetry and thus should be classified according to one of the three symmetry 
classes discussed above. 
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VI. LEVEL REPULSION: CLASSICAL AND QUANTUM ANALOGY 
A. Classical plasma with logarithmic interaction 

One that for the Gaussian invariant ensemble one can rewrite the JPDF Eq. (fTQ| in the following way: 

JP(H)ocexp[-/?£], C = -^ln\E n -E m \+a'J2 E n, (27) 

n>m n 

where we introduces the Dyson symmetry parameter: 

{1, for real symmetric H 

2, for complex Hermitean H (28) 

4, for real — quaternionic Hermitean H 

Note that by a proper choice of energy units the parameter a' = a(3~ 1 can be set equal to \ which will be always 
assumed throughout the lecture notes. Thus there is only one important parameter (3 in the classic WD random 
matrix theory. 

Looking at Eq. (j2"T|) one concludes that the PDF in the (f, E) representation coincides with the partition function 
of classical particles repelling each other logarithmically, in a harmonic confinement potential. The Dyson symmetry 
parameter /3 plays a role if an inverse temperature. 

The above derivation which lead to Eq. (|27[) can be repeated for an arbitrary invariant RME. The corresponding 
energy functional C of the logarithmically repelling particles will differ from Eq. (|27p only by the confinement potential 
which will be no longer harmonic but rather (3~ l V(E n ). The basic property of the PDF which dependends only on 
the set of eigenvalues E n (but not the eigenvector variables f) is retained for all the invariant RME making the 
corresponding eigenfunction statistics trivial. 

This is no longer true once the invariance under basis rotation is broken. The latter circumstance is what makes 
non-invariant ensembles difficult to solve but at the same time having a rich variety of eigenfunction statistics. 



B. Quantum analogy 

Besides the analogy with logarithmically repelling classical particles at finite temperature /3 _1 living in one dimension 
(Id) there is also an important analogy with the system of quantum particles in Id. To facilitate this analogy let us 
remind that the Jacobian in Eq.@ can be expressed as the power of the Vandermond determinant Eq. (|ll[) 

Joz\Af. (29) 

The property of the Vandermond determinant is that 



N d*A N 



E 



n—1 n 



= 0. (30) 



Another property is that it changes sign upon any permutation of two E n and E m . 

These two properties imply that A^r can be considered as the many-body wave function ty({E n }) of the system of 
N free fermions with the Hamiltonian consisting only of kinetic energy: 



N 02 



UK = - (31) 

n=l " 

Moreover, as the energy of the corresponding many-body state is minimal possible for kinetic energy £ = 0, this is a 
ground state of this free fermionic system. 

So we come to the statement that the Jacobian Eq.@ at (3 = 2 is the probability density for the ground state of 
the free fermion system in the entire space. 

J ~ |* (£n)| 2 « H \E n - E^ 13 . (32) 
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Note that the system of fermions in an infinite space is not well defined, as it expands indefinitely. Formally this 
is seen from the fact that the wave function Eg. ([5^|) is not normalizable. In order to fix this pathology one has to 
consider a full probability distribution function Eq.(fT0|) which includes also the confinement potential V(E n ). 
For the harmonic confinement potential the property Eq. (|30[) can be generalized in the following way: 
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(33) 



Now we see that 



(34) 



is an eigenfunction of the free fermions with mass m in a harmonic confinement potential V(E) = ^E 2 . It corresponds 
to a certain positive energy £/v which arises due to confinement of fermions. 

Now suppose that this property is valid also for arbitrary (3 and check that the wave function 



*0,/3({£n}) = II \ E " ~ E m f /2 sgn{E n - E n 



)e ^ 



(35) 



is the eigenfunction of a certain Hamiltonian. Note that the coefficient m can be done arbitrary small by a proper 
choice of E n units. So, for simplicity of further derivation we consider the case m — > 0. 

To this end we take the sum of second derivatives of the wave function applying the kinetic energy operator Eq. (f3Tj) 
to Eq. ([35|) with m — ► 0. The result appears to be proportional to fyo,p({E n })'- 



(36) 



Thus at any f3 ^ 2 the system of fermions equivalent to an invariant random matrix theory is interacting with the 
interaction Hamiltonian: 
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(37) 



Now if we use the property of the Vandermond determinant: 
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we finally obtain the total Hamiltonian of an equivalent system of fermions: 



(38) 



H = - 
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2 V2 -J^ m (E n -E m f 



(39) 



This is the celebrated Calogero-Sutherland Hamiltonian^ with the inverse square interaction. For (3 = 2 the interaction 
constant vanishes and the entire level repulsion is due to fermionic nature of the fictitious particles. For (3=1 there 
is some attraction on top of the free fermionic mutual avoiding, while for [3 = 4 the interaction is repelling. This 
additional interaction explains why the level repulsion for (3 = 4 is stronger then for (3 = 2 and for (3 = 1 it is weaker 
than for (3 = 2. 



VII. PLASMA MODEL AND THE WIGNER SEMI-CIRCLE 



The model of classical particles in one dimension with logarithmic repulsion Eq. (j27]) can be represented by a 
continuous energy functional: 



C=- 



i / dE J dE' q{E) q(E') In \E-E'\ + ^ j dE q{E) V(E), 



(40) 
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expressed through the exact density 

q(E!) =Y,8(E - E n ). (41) 

n 

Now we make two assumptions: 

(i) replace g(E) by an ensemble average value p(E) and 

(ii) neglect the thermal fluctuations by minimizing the energy functional Eq. (|40|) (with g replaced by p) instead of 
computing the partition function 

E e_/3£ ( 42 ) 

As a result one gets a kind of mean field approximation which is justified by the long-range, logarithmic nature of 
interaction. 

Minimizing Eq. (|40|) with respect to p(E) and differentiating both sides with respect to E one obtains: 

+ °V)^=^/(£). (43) 
Hy 1 E - E' (3 dE J K ' v ' 

The physical meaning of this equation is very simple: the force acting upon the given "particle" from all other particles 
should be balanced by the confining force. This is the condition of the plasma equilibrium. 

From the mathematical viewpoint Eq. (|43[) is a strongly singular integral equation. Its solution is well known&. For 
an even function V(E) = V{—E) it reads: 



Po{E ) = -1 1° (44) 

mK ' 7T 2 V J_ D y/D 2 -E' 2 E'-E 1 1 1 



where the principle value of the integral is assumed in Ed.(|43p and Ea.(|44p. namely 

11/1 1 



E'-E 2 \E' - E-iO E' - E + iO 



(45) 



This definition allows to make an analytic continuation of Eq. (|44p for E in the complex plane with the cut along 
the real axis with \E\ > D, where the bandwidth D should be chosen from the condition that the total number of 
eigenvalues is equal to the size of matrix N: 

D 

p (E)dE = N. (46) 

-D 

Namely, p{E) can be represented as a sum of a function p+{E) which is regular in the upper half-plane ^E > and 
a function p-(E) which is regular in the lower half-plane ^sE < 0: 

ME) = M£) + M£), »W - ^ ^^f_ B y& W -f T .u ' (47) 



It is important that along the cut |J5| > D the analytic function \JD 2 — E 2 = ±i\/ E 2 — D 2 has different signs just 
above and just below the cut. This means that for \E\ > D 

p + (E) +p _{E) = -- ^/W^IO 2 [° dE' S(E E') = 0. (48) 

7T J-D VD 2 - E' 2 

On the other hand, for real E beyond the cut (\E\ < D) one obtains: 

p + (E) - p_{E) = ^ JW^W 2 [ D -jffL= S(E -E') = - f(E). (49) 



2tt 2 J_ d y/D 2 - E' 2 v ' it 

Now we are in a position to check that Eq. (f4"4")) is really a solution of Eq. (|4U)) for real E beyond (\E\ < D) the cut. 
Indeed, the integral over the real axis in Eq. (|4"51) can be closed either through the upper complex half-plane of E 1 or 
through the lower half-plane. We use the first option for the part containing p + (E') and the second option for the 
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part containing p_(E'). Each of the two contour integrals allows for the evaluation using the residue theorem. Then 
omitting the terms which do not have poles in the corresponding half-plane we obtain: 

+CX> po(jB ')^_ = i / P±^l dE ' + l f dE' = -ni[p + (E)-p_(E)]=f(E). 

^ m > E _ E , 2 J uppcl E-E' + iO 2 7 lowcr E - E< - zO [P+{> Pyn JK) 

This concludes the proof that Eq. (|44| is indeed a solution of the integral equation Eq. (|43| . The beauty of the proof 
is that it is based only on the analytic properties of the solution. 

For the Gaussian ensemble where f(E') = E 1 the integral in Eq. p^|) is actually independent of E (show this 
using the definition of the principle value of the integral) 

D — — (° 1 dE > _ 

- D VD 2 - W E> -E~ J_ D VL> 2 - E /2 ~~ E ~ n ' ( ' ,U1 
and the average density is the celebrated semi-circle: 

p (E) = - \/2N - E 2 . (51) 



VIII. PROBABILITY OF HAVING A HOLE IN SPECTRUM AND THE WIGNER SURMISE 

One of the most popular statistics of eigenvalues of complex quantum systems is the the level spacing distribution 
P(w): the probability density to have a level at a distance to from a given level and no other levels between them. For 
u> much smaller than the mean level spacing A = p 1 , it is improbable that in between of the two close levels there 
is yet another one or several levels. Then the requirement of having no levels in between of the two is unimportant 
and the leading term in P(to) is the same as in the two-level correlation function R(oj) oc w' 3 at to <C A. However, 
for ijj 3> A the two statistics dramatically differ: R(to) tends to a constant whereas P(u>) is very small due to a small 
probability to have no levels in between of the two levels separated by a large distance. Basically the P{u>) for oj 3> A 
is limited by the probability of having a hole of the size to in the spectrum. Let us find this probability using the 
plasma analogy. 

As for any fluctuation, the probability of having a hole is given by the energy cost SC of this configuration relative 
to the equilibrium one: 

P(u) oc exp(-p6C). (52) 
One can cast the energy difference in the following way: 

A£ = i f dE f dE' 8p{E)8p{E')\n\E-E'\- \ f ' dE f ' dE' p (E) p (E') hx\E - E'\, (53) 

2 Jc JC 2 J- u /2 J-ui/2 

where the integrals in the first term run over the real axis outside the gap region and in the second term they run 
over the gap region; po(E) is the equilibrium density without the gap and Sp(E) = Pu>{E) — pa(E) with Pui(E) being 
the solution of the integral equation Eq. (|43p with the additional condition that there is a gap for \E\ < to/2. 
The solution with the gap can also be constructed and and looks as follows: 

P K ' 7T 2 ^E 2 - (lu/2) 2 Lli VD 2 - E' 2 E' 2 -E 2 1 ' 

It is important that for the steep confinement V(E) ~ \E\ a (a > 1) there is a scale separation, namely the integral in 
Eq. (fS~i)) varies slowly as a function of E with the typical scale of D ~ N 1 ^ . In the large N limit one can disregard 
this dependence and consider 

\E\ 

Puj(E) = po ■ 11 (55) 
^E 2 - (w/2) 2 V ; 

One can immediately recognize the gapped density of states with the square-root divergency near the gap edges similar 
to the one for a BCS superconductor. 
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Now by making a re-scaling E — > sir, E' — > sx' and observing that the double integral in the first term is convergent 
f° r Pu{E) of the form Eq. ([55|) we immediately obtain that AC = lo 2 (a + b\nui). More detailed inspection show that 
the coefficient b = 0. The final result for AC reads: 

AC= ^(po^) 2 «0.62(p ^) 2 . (56) 
Id 

This implies that the spacing distribution function for large level separations s — {lo/ A) 3> 1 is given by: 

P(s)ocexp(-^/). (57) 

Note that the "Gaussian" form of P(s) has nothing to do with the quadratic confinement potential (Gaussian invariant 
ensemble). In fact P(s) has the same asymptotic form Eq. ([57|) for all steep confinement potentials, 
check that it has the same form for the confinement potential V(E) = E 4 . 
Finally we mention a famous interpolation formula for P(s) known as the WignerSurmise: 

P(s)^A(0)s^exp[-B(0) s 2 ], s = ^. (58) 

The coefficients A{0) and B(j3) are found from two conditions: the normalization to the total probability 1 and the 
condition that the mean level spacing in the units of s is one: 



/ P(s)ds = l, / sP(s)ds = l. 
Jo Jo 



(59) 



These conditions result in: 



4 



0.78 (0.62) = 1 



A(0) = 2B(0^ +1 t(^- + i) , B{0) = { | w 1.27 (1.24) = 2 (60) 



2 



64 
9tt 



2.26(2.48) = 4 



For comparison we give in the brackets the exact values of B CK&ct (0) = ir 2 0/16. One can see that they are rather 
close to the approximate values of the Wingner Surmise, especially for = 2. 



IX. LEVEL COMPRESSIBILITY, NORMALIZATION SUM RULE AND NORMALIZATION ANOMALY 

The two-level correlation function (TLCF) is formally defined as a correlation function of the exact density of states 
Eq.flH]): 

Rn(E, E') = ^ly^ly = R-\E)5{E -E') + l + Y N (E, E>). (61) 



The 5- function in Eq. lfBTj) is the self-correlation coming from one and the same level n in the sum in Eq. (f4Tj) . The 
1 term gives the asymptotic value of TLCF at energy separations \E — E'\ ^> A when the average of two densities 
of states can be decoupled. The function Y(E, E') gives then a regular contribution to the TLCF which decreases to 
zero as \E — E'\ increase. 

There is an important normalization sum rule that applies to the TLCF. Indeed, consider 



dE' p{E) p(E') R N (E, E') = J dE 1 (g(E) g(E')) = (g(E) J dE 1 g(E')) (62) 
The total number of states is equal to the number of degrees of freedom N and does not fluctuate: 

J dE' g(E r ) = N. (63) 



This normalization condition leads to 



J dE'p{E)p(E')R N {E,E') = N(g{E))=Np{E), 
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which implies that 

/+oo r + cc 

dE' p(E') Y N (E, E') = -1, / dE' p(E') [R N (E, E') - 1] = 0. (64) 

-oo J — oo 

This is the normalization sum rule. 

Note that the sum rule can only be proven if N is finite and integration in Eq. (l64j) are extended over all energies. 
Taking the limit N — > oo could be a dangerous procedure as in this case one has to worry about the commutativity of 
the limits N — > oo and limits of integration— ■> oo. The sum rule is certainly satisfied if the limit N — > oo is done after 
doing the integral. A simple example below shows that it can be violated if the limit N — > is taken prior of doing the 
integral. 

Consider an ensemble of diagonal random matrices with independently fluctuating random elements each having a 
distribution 

Pie ) - ( N ~'' I 6 "' < N l 2 (65) 
ne " j -\ \e n \>N/2 [bb) 



The TLCF for this ensemble can be computed straightforwardly: 

/N/2 7 pN/2 j pN/2 j 

*S / zE^L S(E - e n ) 5(E< e m ) + W ^S(E- e n ) S(E> - E „), 
„..,„ -N/2 A J_ N/2 N - J_ N/2 N 



so that for \E\ < N/2 and \E'\ < N/2 one obtains: 



N 2 - N 1 
R(E,E')=S(E-E')+ N2 , p(E) = l, Y N (E,E') = --. (66) 



At a finite N the normalization sum rule Eq. ([64|) is obviously fulfilled: 



N/2 , i 

dE' --)=-!. (67) 



N/2 



N 

However, if one takes the limit N — > oo in Eq. (|66|) before integrating, one obtains Yqo = 0, and the normalization sum 
rule will we violated: 



f N/2 

/ dE' lim Y N (E,E') = 0. (68) 

J-N/2 N~*oo 



This mechanism of violation of sum rules in the thermodynamic limit is called the anomaly and is well known in the 
field theory. 

A remarkable property of the Wigner-Dyson level statistics is that in this case (as well for all invariant RM 
ensembles, even with shallow confinement potentials) the normalization sum rule is not violated and the anomaly 
does not occur. Let us show how this property follows from the plasma model Eq. (|40p . To this end we note that 
the density-density correlation function is given by the variational derivative of the mean density with respect to the 
confinement potential: 

(g(E) g(E')) - (g(E)) (g(E')) = « pi [R^E E') 1], (69) 

where we assume that the mean density p[E) ~ po does not change much at a scale of the mean level spacing A. In 
this case one can approximately consider the TLCF as a function of the energy difference. Then integrating by parts 
in Eq. (|44[) and neglecting the energy dependence of the square roots we obtain: 

S P (E) _ 1 1 



SV(E') tt 2 /3 (E — E') 2 ' 
where the regularization 

(E - E'y 2 — ► i [(E — E' + iO)- 2 + (E-E'- iO)- 2 } 
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is assumed. Using this regularization one can immediately check that 

/+oo 
[R oo (E)-l]dE = 0, 
-oo 

as it is required by the sum rule Eq. ([64|) at p(E) = const. Thus we see that for the plasma model doing the limit 
N — > oo and doing the integral commute. 

The absence of the anomaly is related with the incompressible character of the system of logarithmically repelling 
particles. This is the reason why approximations made in deriving the plasma model Eq. (|40|) did not affect the regular 
fulfillment of the normalization sum rule. Below we define the level compressibility and show that it is zero if the 
normalization sum rule is not violated. To this end we define the level number variance: 

E(n) = (n 2 ) -n 2 , (71) 

where n is the fluctuating number of levels in an energy interval Se that contains on the average n levels. Writing 

rE +6E/2 
I Eq-SE/2 



pE„+SE/2 p 

/ p(E)dE= / p{E)dE 

JEn-SE/2 JSk 



we obtain: 



S(n) = j p(E) dE I p(E') [R N {E, E') - 1] dE'. 



Now we assume that N is large, the confinement is steep and thus D ^> \5E\ ^> A. This allows to consider p(E) = po 
and Rn(E, E f ) » Roo(E — E'). Then one can integrate over E at a fixed E — E' and arrive at: 

/n rfi 
[Roo(s)-l]ds- \s\[R oc (s)-l}ds, (72) 
-n J — n 

where s = (E — E')po. The level compressibility is defined as 

dY P" 

X(n) = ^= / [Roo(s)-l]ds. (73) 
an J_ n 

In the case where the anomaly does not occur we have 



[Roo(s) - 1] ds = lim / [R N (s) - 1] ds, 



which in view of the normalization sum rule Eq. (|64j) implies incompressible nature of the system of energy levels: 

lim X {n) = 0. (74) 



The opposite is also true: if Eq. ([74|) is fulfilled, the limits n — > oo and N — >■ oo commute and there is no anomaly. 

Now we use the absence of the anomaly and the normalization sum rule to compute the level number variance at 
n 3> 1. To this end we cast Eq.([72|) in the following way 

S(rl) = n 



[R 00 (s)-l]ds-2n [R 00 (s)-l}ds-2 \s\ [^(s) - 1] ds, 

-oo J n JO 



The first integral vanishes because of the sum rule Eq. ([M|) . the second term requires only the knowledge of TLCF at 
large distances s ^> 1 and can be computed using Eg. ([70)) and appears to be a constant of order 1. The third integral 
is logarithmic and this allows to compute the leading logarithmic term also using Eq. (|7D|) which is valid at s > 1. 
Cutting the logarithmic divergency at s ~ 1 we obtain for n> 1: 



S(n) = 4tt Inn + 0(1). (75) 



This variance is considerably smaller than for independently fluctuating levels (diagonal RME) where it is distributed 
according to Poisson law: 

T,(n) = ft. (76) 
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X. ORTHOGONAL POLYNOMIALS AND ENERGY LEVEL STATISTICS FOR /3 = 2. 



As has been already mentioned level statistics in the Gaussian invariant RME with (3 — 2 can be exactly mapped 
onto the system of non-interacting fermions in one dimension in the parabolic confinement potential V(E) = E 2 . 
Let us check this statement. To this end we recall that the ground state many-body wavefunction of non-interacting 
fermions is the Slatter determinant: 



*({En}) = 



fo(Ei) 
<Pi(Ei) 



tfo(E 2 ) 
<fii(E 2 ) 



<Po(E N ) 
(pi(E N ) 



ip N -x{Ex) ip N ^i(E 2 ) ... tp N ~i(E N ) 



(77) 



The one-particle eigenf unctions f n (E) = H n (E) e E I 2 in the parabolic confinement potential V{E) obeying the 
Schroedinger equation: 



\ ^2 ^(E) + l -E 2 <p n (E) = £ n <p n (E), 



(78) 



are related to the Hcrmitc orthogonal polynomials H n {E). These are the polynomials of the n-th order satisfying the 
orthogonality relation: 



exp[-V(x)] H n (x) H m (x) = h n S n 



V(x) 



1. 



It is important that the Hcrmitc polynomials obey the three-term recursive relation: 

H n+ i(x) = 2xH n (x) - 2nH n -i(x). 
Using this relation one can show that the Slatter determinant Eg. ([77]) reduces to: 



<£o(-Ei) 


<Po(E 2 ) . 


■ Vo(En) 




1 


1 


1 




r 1 N i 


<Pi (Ex) 


tpi(E 2 ) 


ipi(E N ) 


= const 


Ex 


E 2 . 


En 


exp 


ip N -i(Ei) 


ip N -i(E 2 ) ■ 


■ ¥n-i{En) 




E?- 1 


E f-i : 


■ N 




n=l 



(79) 



(80) 



(81) 



Indeed, the exponential factors e~^ E in all the <p n (E) can be taken out of the determinant using the rule of multipli- 
cation of determinant by a factor which is equivalent to multiplication of all the elements in a column by this factor. 
So we obtain the exponential factor in the r.h.s. of Eq. (|5Tj) . 

Next choosing Hq(E) = 1, H\(E) = x one can find all the other polynomials using the recursion relation Eg. ([50)1 . 
In particular, H 2 — 2x 2 — 2. This polynomial should be plugged in the third line of the determinant in the l.h.s. of 
Eg. ([81]) . Note, however, that the constant term —2 can be omitted as its inclusion is equivalent to an addition of the 
first line to the third line in the determinant which according to the basic property of a determinant does not change 
its value. This process can be continued. For instance in H3 — 2xH 2 — AH\ which stands in the fourth line of the 
determinant one can omit —AH\, as Hi stands in the second line. Then expressing H 2 = 2xHi — Hq = 2x 2 — 1 one 
can put in the fourth line 2x(2x 2 — 1) instead of H3. Finally, observing that the term —2x can be considered as a 
linear combination of Ax and the first line in the determinant and omitting this term we conclude that instead of H3 
on can put in the fourth line of the determinant just one term Ax 3 . 

The determinant in the r.h.s. is the famous Vandermond determinant which is equal to: 



(82) 



1 


1 


1 




Ei 


E 2 


En 


— Y\_ (^n — E m ) 


E N-l 


E?- 1 '. 


■ E N 





Now we see that: 



I* ({E n }) | 2 = const exp 



\{{En 



En 



(83) 



This is exactly the probability distribution functions for the eigenvalues of the Gaussian RME with j3 = 2. 
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It turns out that the theory of orthogonal polynomials^ is the powerful method to solve any orthogonal random 
matrix ensemble. What one has to do for that is to generate a set of orthogonal polynomials p n {x) obeying the 
orthogonality relation Eq. (|79| and to be able to compute the large- N asymptotic behavior of the " wavefunctions" : 



<p n (E) = y n {E) = h-WpnWe-vW* 



(84) 



The generation of orthogonal polynomials is possible for any confinement potential V(x) using the Gram-Schmidt 
orthogonalization procedure. According to this procedure one computes the Gram-Schmidt determinant: 



(85) 





Of) 


ai . 


■ a n 


G n — 


ai 


a 2 . 


■ a n + l 




a„ 


Ofi+l • 


■ 0,2n 



where a n are the moments: 

/+oo 
exp[— V{x)\ x n dx. 
-oo 

Then the orthogonal polynomial of n-th power p n is given by: 

Vn{x) 



(86) 





do 


ai . 


a„ 




Ol 


a 2 ■ 


■ On+l 








G n -i 








O-n-l 


a n . 


• 02n— 1 




1 


x . 


. X n 



(87) 



where q n is the coefficient in front of x n in p n . 

It follows from this generic procedure that any set of orthogonal polynomials p n (%) should obey the three-term 
recursion relation similar to Eq. (|80[) . For the confinement potential V(x) being an even function of x and the choice 
q n = 1 (for the Hermite polynomials the standard definition corresponds to q n = 2") it reads: 



p„+i{x) = xp„(x) - C„+i p n _i(a;), 



G n G n —2 , , , , 

Gn+l = > POW = 1 1 PlK X ) = X - 



We note that this recursion relation generates orthogonal but not ortho-normal polynomials. The price of having 
q n = 1 is that the normalization constant h n in Eq. (|79j) is not unity and is related to the coefficient C n as follows: 



n+l 



hn — ho JJ C m , ho — 



(89) 



m=2 



The recursive relation Eq. (|88|) appears to be the most convenient way of generating orthogonal polynomials for 
any confinement potential. It works also for non-classical polynomials for which there are no second-order differential 
equations (similar to the Schroedinger equation Eq. (|78p in the case of Hermite polynomials) which the " wavefunctions" 
Eq.([8~4)l should obey. 

The efficiency of the orthogonal polynomials in the problem of level statistics is largely due to the Christoffel-Darboux 
formula: 



N-l 



K N (x,y) = 22 (p n (x)(p n (y) 



h N y>jv_i(x) <Pn{v) - <fiN-i(y) <Pn(x) 



flN-l 



y-x 



and a relation Eq. 



(90) 
between 



This formula can be proven by induction using the three term recursive relation Eq. 
h n and C n . 

The Christoffel-Darboux formula is important because the mean density of states and the two level correlation 
function at x ^ y are given by: 



N-l 



p(E)=Y,v n (x) 2 = K N (E,E), 



(91) 



71=0 
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JV-l N-l 



p(E)p(E')R N (E,E')=J2 E H(E)^ m (E')-ip n (x) Vn (y)^ m (y)^ m (x)] = K N (E,E) K N (E',E') - K%(E,E'). 



n—0 m— 

(92) 

Eqs. (|9ip . (|92p can be proven formally without any reference to systems of non-interacting fermions. However, it is 
instructive to see how they follow from the fermionic second quantization formalism. Indeed, the density of non- 
interacting fermions and the density-density correlation function are given by: 

(0|*t(aOfc(aO|0), (0|fct(a:)*(a:) &(y)*(y)\0), 

where (0|...|0) is the quantum-mechanical averaging over the ground state. According to the rules of second quanti- 
zation the fermionic field operator ty(x) is given by the expansion over the single-particle wavefunctions: 

N-l 



n=0 



where c„ and cj, are the fermionic creation and annihilation operators obeying the anti-commutation relation 



The averages over the ground state can be computed using the Wick theorem: 



]clc m \0) = 5„ 



C n 1 c n 2 C n 3 c n i 



|0) = S nin2 S n 



(93) 



The two terms in Eq. (f92|) follow from the two terms in Eq. (l93|) which correspond to two possible parings of c? and c. 
Eq. (l92p can be conveniently represented in terms of the determinant: 



p(E)p(E')R N (E,E r ) = 



K N (E,E) K N (E,E') 
K N {E',E) K N {E',E'] 



(94) 



Using the fermionic analogy and the Wick theorem one can prove that any multi-point level density correlation 
function can be represented in the form of a similar determinant: 



p{E 1 )p{E 2 )...p{E n )R{E 1 ,E 2 ...E n 



Kn{E\, Ei) K N {E 1 ,E 2 ) ... K N (Ex,E n ) 
K N (E 2 ,E 1 ) K N {E 2 ,E 2 ) ... K N (E 2 ,E n ) 

K N {E n ,E x ) K N {E n ,E 2 ) ... K N (E n ,E n ) 



(95) 



We see that the analogy with non-interacting fermions allows to express any multi-point level density correlation 
function in terms of only one single kernel K]y(x,y). The latter according to the Christoffcl-Darboux theorem is a 
product of only two "wave functions" which require the knowledge of only two orthogonal polynomials pn(x) and 
Pn-i(x). Thus the problem of energy level statistics is reduced to the problem of finding the asymptotic behavior of 
orthogonal polynomials of high order. 



XI. WKB QUASI-CLASSICAL APPROXIMATION AND THE ONE-DIMENSIONAL "WIGNER 

CRYSTAL". 

The semicircle law for the mean level density follows immediately and trivially from the free-fermion representation. 
Indeed, the density of one-dimensional fermions is directly related with the Fermi-momentum pp by: 

ff=P> ( 96 ) 

When the density varies slowly at a scale of the Fermi wavelength, one can apply Eq. (f9T>|) locally thus relating the 
local Fermi- momentum Pf(x) with the local density p{x). The local Fermi momentum corresponds to the momentum 
of the highest occupied state in a parabolic potential 



PF = Pn{ x ) — \J 2™£'kin = V ^Sn — x 2 = \J 2N + 1 — x 2 , m = 1. 



(97) 



15 



Then the local density is obtained immediately from Eq. (f96j) : 

p(x) = -Pn(x) = - V2N + l~x 2 . (98) 

7T 7T 

This is the celebrated semicircle law. 

In order to obtain the density-density correlation function, or the two-level correlation function one should work a 
little bit harder determining the large N asymptotic behavior of tp n and applying Eq. (|92| . 

In the case of Hermite polynomials the problem of the large- N asymptotic behavior can be solved quite easily. 
The reason is that there is the second-order differential equation (the Schroedinger equation) which the wavefunctions 
<Pn{x) must obey. As is well known the solutions to the Schroediger equation corresponding to large quantum numbers 
bear the properties of classical motion in the corresponding potential. In quantum mechanics this corresponds to the 
" quasi-classical" , or WKB approximation^. 

According to this approximation the wave function at x 2 < 2N is proportional to: 

{cos [f£ pn(x') dx'] N is even 
(99) 
sin [ Jq Pn(x') dx'\ N is odd 

At large TV we obtain: 

p N (x)^V2N, [ p N {x')dx' « x V2N. 
Jo 

Then the kernel K^(x,y) is easily calculated using the Christoffel-Darboux formula Eq. (|5D|) : 

sin (V2N (x - y)\ 

Kn(x, y) = const — , Kjyfx, x) = const \2N. 

x-y 

The normalization constant const = 1/tt is most easily found from the comparison of Kn(x, x) — p(x) and the semi- 
circle mean level density Eq. (f5Tj) which at large N reduces to p(x) « \/2N/ir — pa. Now, introducing mean level 
spacing A = p^ 1 w tt/\ / 2N we arrive at: 

K N (x,y) -f K(s) = po S -^^-, s=^-, 7V->oo. (100) 



7T.S 

The two level correlation function Eq. (j!)2")) is then equal to: 



A 



r, / \ x sin 2 (7rs) , 1 cos(27rs) 

Roo (x, y) = S(s) + 1 - = l - _ + ^L_i. (101) 

One can see that the TLCF given by Eq. (|10ip has all the asymptotic limits right. It is proportional to s 2 at s <C 1 
and its envelope corresponds to Eq.ffOl) for s ^> 1. In addition to that it obeys the normalization sum rule Eq . ([64| . 
However, in Eq. p01[) there is a term that oscillates with the period of the mean level spacing A. This term evades 
the consideration based on the 2x2 matrix (small s<l) and the effective continuous plasma model (large s ^S> 1). 

Let us discuss the physical meaning of this term using the plasma model analogy but without the continuous 
approximation. It is well known that plasma of particles with the long-range repulsion in a confinement potential 
tend to develop a crystal order known as Wigner crystal. Such Wigner crystal of electrons have been observed on 
top of the helium surface. Our case is special, as it is one-dimensional. According to the Mermin theorem the crystal 
order cannot survive in one dimensions at a finite temperature because thermal fluctuations destroy the long-range 
order. However, local crystal order may exist. The last oscillating them in Eq. ljlOip reflects exactly this order. The 
short-range nature of this order manifests itself in the fast s~ 2 decay of oscillations at large distances. 

So far in this section we have considered the (5 = 2 case. As in the plasma analogy (3 plays a role of inverse 
temperature, one would expect the oscillating term to decay slower for (3 = 4 and faster for (3 = 1. This expectation 
is in fact true. 

One can show using the more sophisticated application of the orthogonal polynomial machinery that in the limit 
N — > oo the two- level correlation functions for the orthogonal ((3=1) ensemble i?^ th and that for the symplectic 
((3 = 4) ensemble i?^ mp can also be expressed in terms of the kernel K^r(s), Eq. l|100[) : 

Y^(s) = -K\s) - d -^- ^ K(x) dx, (102) 



1G 



l™( S /2) = -X 2 ( S ) + ^M / KU),i,. (UK!) 



o 



The s/2 in the argument of T^ mp (s/2) appears because of the Kramers degeneracy: for the same total number of 
levels N the mean level spacing between doubly degenerate levels is two times longer. Accordingly, the S(s) function 
in R(s) enters with the pre-factor of 2. The asymptotic behavior of these functions for s ^> 1 is the following: 



^) = -i + ^, (104) 



ysymp r x_ 1 C0S(27T S ) cos(47T s) 

One can see that the leading oscillating term decrease as s~ 4 /^ as was expected. In addition to that, in the symplectic 
ensemble = 4, the sub-leading second-harmonic term appears which was absent to all orders in 1/s for (3 = 1,2. 

The spectral correlations of a quantum system show up in the timc-dcpcndcncc of response to external time- 
dependent perturbations. For such applications one need to know the Fourier-transform F(t) of the two-level corre- 
lation (cluster) function Y oc (s). It appears to be amazingly simple for the unitary ensemble f3 = 2: 



t\-l, \t\<l 
0, \t\ > 1. 

with the jump of the first derivative at t = 1 that leads to the oscillations with the period 1 which amplitude decreases 
as 1/s 2 . For the orthogonal ensemble [3 = 1 there is a jump only in the third derivative: 

f2|t|- l-li| ln(l + 2|*|), |*| <1 
portho J / , , \ n n7\ 



7 ln(^j, |*| >1. 

For the symplectic ensemble (3 = 4 there are two singular points: |*| = 1 and |*| = 2 which correspond to two oscillating 
terms in Eq. (|105[) : 

F s ymp = R 1*1-1- \\t\ In |1 - |*||, |*| < 2 (1Q8) 



XII. WIGNER-DYSON LEVEL STATISTICS AND THE LUTTINGER LIQUID. 

The large- s asymptotics of the two-level correlation function containing both the non-oscillating and the oscillating 
terms which decay as a certain power-law can be written in a compact form which involves only one single function 
G(s): 

Y ^( s ) = -4^2 + cos ( 2 ™) ^ (1° 9 ) 

Y£ tho ( S ) =~-^^^ + 2cos(2^)e 2G ( s ), (110) 

^ ymp (*) = -^2 ^r- + ^ cos ( 2 -) eG(s)/2 + 1 cos ( 4 -) e2G(s) - (in) 

where 

G= -ln(2^ 2 s 2 ). (112) 

It turns out that the function G(x) is proportional to the equal-time correlation function of a free bosonic field ^(x, r) 
in the too dimensional space- (imaginary)time, which arises as a bosonized version of the Calogero-Sutherland model 
Eq. (|39[) of interacted fermions. More generally, a great number of models of interacted electrons in one dimension 
fall into the universality class of Luttinger liquid^ which is characterized by a certain correlation functions at large 
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separations in space and/or in time. All of them follow from the fact that the fermionic operator t) can be 
represented as 

V(x,t) = Re ikpx + Le~ ikFX , R, L = -= exp[±i$ fl L (x, r)], (113) 

V7T 



where $(x, t) = 5 [&r(x, t) + $>l(x, t)] is the free boson field with the action: 

, rl/T p+OO j 

m = 2^K dT dx [( ^ $)2 + {dT ^ )2] = ~jtKT ^ lAk ^ |2 {UJ " + fc2) ' 



(114) 



where 



*(a:,r)= ]T {4 fe)W e < ( w » T+te > + c.c} , (115) 

fe>0 : w„=2irTn 

w n = 2irTn (n are aZZ integers) and k = (2ir/L)m > to avoid double-counting. The action Eq. (|114[) corresponds to: 

,1 . ,9, ttKT 1 

<^l> = — FT< (116) 

It is remarkable that interaction of fermions is encoded in only one single parameter A which is A < 1 for repulsion, 
A' = 1 for the non-interacting fermions and A > 1 for attraction. In other words, with respect to long-distance 
properties the system of interacted fermions in one dimensions (1 + 1 space-time) is equivalent to a system of free 
bosons. The physical meaning of this result is that for systems of the Luttinger-liquid universality class all the 
multitude of effects of electron interaction reduces to dynamics and thermodynamics of the plasmon collective modes. 
The density operator in this representation is given by: 

p(x, t) - po = — d x $(x, t) + Ax cos[2fc F x + 2$(x, t)\ + A 2 cos[4fc F x + 4$(x, t)} + (117) 

where Ak arc the structural constants which are determined from the details of the system at small distances. 

The first term in Eq. (|117p comes from the combination [R + R + L + L] in ^t$. It is analogous to the 8p oc Vu 
term in hydrodynamics, where u(x) is the mass displacement at a point x. The correct evaluation of this contribution 
requires the regularization R + R = R + (x + a, t)R(x, t) (and the similar regularization for L), where a — 1 is the lattice 
constant which corresponds to kp = n. Oscillating terms proportional to e ,' il2lk i p ^ ar i se from the cross combinations 
R + L + L + R in 'J'Tijr, Note that for interacting fermion system there is a vertex correction that involves momentum 
transfer far away from the Fermi points ±kp. As the result, the density operator expressed in terms of the field 
^(x, t) (which contains only momenta close to the Fermi points) is not simply equal to p — but may have higher 
order terms as well: 

p(x, t) = + ci (* t *) 2 + ... (118) 

It is these higher order terms that generate combinations containing higher harmonics like e ±4lfcF2 \ 

Using Eq. (|117[) and Eq. (|114[) one can express the density-density correlation function through the free bosonic 
correlation function. To this end we use the identity valid for any Gaussian field theory: 

( e <n*(«,T) e -m*(0,0)) = exp /_» T ) _ $( , Q)) 



s 

Now observe that for s>l and T — > we have in the thermodynamic limit L — > 00: 

($( s , 0)*(0, 0)) s - (m 0)$(0, 0)) s = - v / d 1 — w T G ( fi )' 

2 Jo q 4 

where G(s) is given by Eg. (jl 12[> . Finally, for oscillating part of the density-density correlator we obtain: 

( e <n*(»,T) e -in*(o,o)) = exp flKn 2 G(s)\. (119) 
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The non-oscillating part is expressed through the second derivative of the correlation function of free boson field: 

^ (d s 4>( s , 0)d s ,$(s', 0)) s = '^ds G(a). (120) 

Eqs. (|119| . ()120p show that the phenomenology of the Luttinger liquid allows to relate the coefficient in front of the 
non-oscillating part of the density-density correlator with the coefficients in front of G(s) in the exponent determining 
the amplitude of the oscillating terms. First of all we fix the interaction parameter K from the amplitude l/(27r 2 /3) 
of the non-oscillating part. Eq. (|120[) suggests that: 

K = \ (121) 

Then a comparison of Eq. (|119p with the oscillating terms in Eqs. (|109|) . (|110[) . (|llip shows that all coefficients k in 
the amplitude of n-th harmonic e KG ^ cos(2itnx) are equal to k = n 2 K with K given by Eq. (|12ip . This is exactly 
the parameter K that corresponds to the Calogero-Sutherland model Eq. (|3"5)) . Thus we have demonstrated that the 
Wigner-Dyson level statistics at large level separations corresponds to the particle statistics of the Calogero-Sutherland 
model at zero temperature. 

One can ask a question: how the Wigner-Dyson ensemble should be deformed in order to retain this analogy with the 
Calogero-Sutherland model also for finite temperatures T ^ 0. The answer is^ that the proper deformation is given 
by the Gaussian non- invariant ensemble Eq.([5]) at large values of the parameter B. The corresponding temperature 
of the Calogero-Sutherland model is^: 

T =m (122) 

XIII. FIELD THEORIES FOR RANDOM MATRIX ENSEMBLES 

In this section we derive the field theory for an arbitrary Gaussian random matrix ensemble. This formalism, known 
as nonlinear super-symmetric sigma-model^ has been first applied to the Wigner-Dyson random matrix ensemble. 
However, its real strength is in the possibility of extension to the non-invariant random matrix ensembles as well as 
to real disordered conductors with diffusive dynamics of particles. 

We start by writing down the expression for retarded (G R ) or advanced G A Green's function in terms of the 
functional integral over complex variables tp n : 

Gnm = ([E± -H]" 1 )^ = ^ J V^pVip* <p n tp* m exp [iS±[<p]] , (123) 

where 

S ± [<p\ = ± y* [E ± 5ij - Hij\ w (124) 

and E± = E ± (ui/2 + iO). The sign ± stands for G R or G A , respectively. There would be no problem to average 
Eq. (|123p over the Gaussian random entries H nm if not the normalization constant (partition function) Z: 

Z = JvtpVtp*exp[iS ± [ip}}. (125) 

which also depends on H nm . With the Z present one has a problem, the problem of denominator. 

There arc different ways of overcoming this problem, e.g. the replica trick. However, here we use another trick, the 
super- symmetry method^. In the core of this method is the calculus of anti-commuting (Grassmann) variables /i m : 

Mm = "Mm Mn, Mn = 0- (126) 

One can define the Grassmann integral with the shortest table of integrals ever: 

Mn dfi n = -j=, J dfin = 0. (127) 
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Since any function /(/in) = f(0) + f'(0) fi n , the Grassmann integration is essentially a differentiation. Now let us 
compute the integral 



Y\ duldni exp 



^ ' a n 



J n d ^ d ^ (-«t) = y n mi^* m/^ ^ = n ~~ 



(128) 



To accomplish this we expand the exponential function to leave only the term that contains a complete set (i r , 
of Grassmann variables and apply the table of integration Eq. (|127p . 

The corresponding integral over the usual complex variables would give the following result: 



JJd^dipi exp 



^2 <Pn °" fr. 



n 



(129) 



We see a remarkable property: the product of the two integrals is equal to 1. This property remains true for any 
Gaussian integrals of commuting and anti-commuting variables. In particular, 



Z^ 1 = J Vfi*Vfiexp[iS ± [fi}}. 
Now the Green's functions can be represented without the denominator: 

ip n ¥4 exp [iS±[ 



where 



S± M = S ± [<p] + S ± [fi]=±Y, 4 [E±S l3 - Ha] ^ 
Here we introduced the super-vectors ip and ip^: 



and the super-measure: 



T>i/j = Dip* "DtpVp* D/x. 



(130) 



(131) 



(132) 



(133) 



(134) 



The action Eq. (|132p and the integration measure Eq. (|134p are super- symmetric, i.e. the commuting and anti- 
commuting variables enter in a fully symmetric way. The super-symmetry is however broken in the pre- exponent 
in Eq. (|131[ ). as it depends only on the commuting variables. 

Now when the problem of denominator is solved by the supersymmetry trick, the next step is to average over the 
Gaussian ensemble of Hy. To this end we write: 



Averaging of the r.h.s. is done independently for each pair of i,j using the identity: 



r.h.s.--^-\H ij \ 2 



1 



From now on for simplicity we will consider the case (5 = 2. Then 



Ha ± 4^ ) [ h ti ± 1 -a 13 4^ ) - ^ 4^H 



dHijdH*^ exp [ r.h.s. — |iJ 



A, 



J dU, dU: s exp (--L iHij^ exp (-^L 4^4^ , 



where 



H i:j = Hij ± -ipjipj. 
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The simplicity of the case (3 = 2 is that H^ belongs to the same manifold of complex numbers as H^ , so that one may 
replace in the integral over the entire manifold Hij — > Hij. Thus on the right hand side we obtain the normalization 
integral for the random matrix ensemble averaging. So we obtain for the disorder average: 



cxp l^fi ^jHi^j ^ = exp ^-i A v ^Hj^i^j 



(135) 



Now we define the super-matrix 



and the super-trace: 



*®w=te » - a z) d36) 



- ("',1 %■ ) ^ (l:!7. 



Then Eq. (|135p can be conveniently rewritten as 



Finally, the averaged Green's function can be represented as follows: 
(G*L A )=TiJv^ Vn <p^exp[-F[Q]], F[Q] = T^i^STr^ + i^^- STr[QiQj], Qi = ^®^(138) 

Thus we derived the deterministic field theory which is equivalent to the Gaussian random matrix ensemble and is 
suitable to compute the average Green's function. One can see that the matrix of variances 

A ij = (\H ij \ 2 ) 

plays a role of the coupling constant ("coupling matrix") in the corresponding action F[Q]. 

The filed-theory representation Eq. (|138[) can be extended to consider the averaged product of (G R G A ) which is 
necessary to be able to compute the two-point correlation functions. To this end, one introduce the double set of 
commuting and anti-commuting variables: one for G R and another for G A and then repeats with minor modifications 
all the above steps: 

(G*n (E + ~) G A m , (i5-|))= <rf ^ <p A , <p% exp [-F[Q]] , (139) 

where the action S[Q] takes the form: 

F[Q] = -iEY, STr[Q t ] - t ^±^. ^ STr[AQj] + \Y, A » STr KM> Qi = ^i® (140) 



with 



f<p R \ /l 



* = Up R * n R * ~ip A * -li a *) , * 



,010 . 

A = A, 1. I o o _! o I H41I 

v0 -1, 



There is another, in some sense dual, field-theory representation similar to Eq. (|139[) . In contrast to Eq. (|139[) it involves 
the inverse coupling matrix (A )y rather than the variance matrix Ay. In order to obtain this representation one 
makes the Hubbard- Stratonovich transformation: 

eatpj—^^STlrK^i® »*,•)] | = JvPcxp J - ^ Ar. 1 STr^P,] + i £ STt[(9 t ® #;)P;] j , (142) 
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where Pi is a super-matrix field. 

If one substitutes the last term in Eq. p40j) for Eq. (|142p the remaining integral over the \& fields is Gaussian which 
can be done using a generalization of Eqs. (|128ll - (1129|) for the case where the super-matrix Ki is not proportional to 
the unity matrix (for which STr[lnK;] = 0): 



JtW exp^iKiVj} =exp|-^STr[ln^]| 



(143) 



Now we are in a position to write down the full action of the dual representation: 

F[P] = J2 A il Vto[PiPj] + E STr[ln(S - P % + ( w /2) A)]. (144) 

ij i 

In order to appreciate the possibilities this representation is offering and also for further simplifications of Eq. (|144|) 
we compute the inverse coupling matrix A^ 1 for the important case of the banded random matrix ensembles where 
the variance matrix A,-j is given by Eq.(|U). 

The matrix element of a matrix inverse with respect to a matrix Ay = A(i — j) is given by: 



.1,/= ±G(k)e^\ G{k) = 



-, -1 



E A{k + 2nm) 



m— — oo 



r + oo 

A(k) = I drA(r)e~ tkr . (145) 



Note that the summation over the reciprocal lattice vector 2-nm is important. 
For the case of the exponential A(r) = e~' r '/ B , we have: 

- 2B V 2B Sinh( ^ ) 

{ > 1 + S 2 fc 2 ' ^ l + B 2 (fc + 27rm) 2 cosh(i) -cosfc' 

m— — oo \ rs j 



One can see that for large B ^> 1 the inverse coupling matrix A- 1 is extremely simple: 

G{k) » B (1 - cos k), A- 1 » | (2% - «y w+1 - 8ij-i), Ao" 1 = ^ 4^ » ^- (146) 

i 

Eq. (|146j) shows that locally it is just the lattice second derivative. However, the pre- factor in front of it is large ~ B 
and it becomes infinite for the Wigner-Dyson ensemble. This means a large cost of variation of Pi in the space. Let 
us consider Pi = Pq being independent of i to the first approximation. Then the action becomes: 

iV" 1 F[P ] = 1 STrP 2 + STrln(£ - P ). (147) 



Here we also neglected a term proportional to uj which is legitimate as long as uj <C yB. Minimizing the action with 
respect to Pq we obtain a saddle-point equation: 

P (E-P )=Ao/2. (148) 
The solution to this saddle-point equation is degenerate: 

Pv = \(E + iQ v/2Aq - £ 2 ) , (149) 

where Q is a super-matrix obeying the constraint 

Q 2 = 1. (150) 

Another constraint comes from the requirement of super-symmetry: 

STrQ = 0. (151) 

Indeed, a super-matrix Q obeying the constraint Eq. (|150p after the diagonalization must contain only ±1 on the 
diagonal. The super-symmetry (the symmetry between commuting and anti-commuting variables) implies that it 
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must have the same diagonal elements corresponding to fermionic and bosonic variables of the given type (R or A). 
But then necessarily STrQ = 0. 

Now take into account (as the first order expansion in uS) the term in Eq. (| 1441) proportional to the energy difference 
uj and allow for slow variations of the field Q = Qi in space which do not violate the saddle-point condition Eq. (|150p . 
Then we obtain plugging Eq. (|149|) into Eq. (|144p : 

F[Q] = ~ (irpAo) 2 J^Ar. 1 STriQiQj] - i^f £ STr[AQ 4 ], (152) 

ij ~ i 

where 



P = Po(E) = ^ I -^, A (153) 



This is the celebrated action of the nonlinear a mode&^. 

For the particular case of banded random matrices with Aij given by Eq.(j4]) we obtain: 



F[Q] = -Dj2 STr K^ - ^+i) 2 l - ^ E STr [ A( 2i]> ( 154 ) 

i i 

where D = \B (j3 — ^-J ~ B 2 . The continuous limit of this model is the diffusive nonlinear er-model: 

F[Q] = -D J dx STr[(VQ) 2 ] - i^- J dx STr[AQ(x)], (155) 

which was originally derived by Efetov^i to describe the crossover from the diffusive dynamics to the Anderson 
localization in the quasi-one dimensional multi-channel disordered wire. This demonstrates the isomorphism of the 
problem of quasi-one dimensional localization and the problem of banded random matrices^. 

Closing this chapter we note that the derivation of Eq. (|152|l from the exact Eq. (|144ll requires the saddle-point 
approximation Eq. (|148[) . Thus Eq. (|152[) is justified only if the energy cost of space variations of Qi is high. This 
happens when the variance matrix Aij has a form of a banded matrix which is approximately constant at \i — j\ < B 1 
where B ^> 1. For the Wigner-Dyson ensemble Aij = 1 and the bandwidth is maximum possible B ~ N (in 
particular, Aq = N). In the limit N — > oo all spatially varying configurations of the field Q (non-zero modes) are 
strictly forbidden. Neglecting them we obtain the zero-mode nonlinear sigma-model-- which describes the statistics 
of energy levels in fully chaotic quantum systems of confined geometry (quantum dots): 

F WD [Q] = -i^-STr[AQ], s = ^. (156) 



XIV. HOW TO COMPUTE OBSERVABLES: SEMI-CIRCLE LAW FROM THE SOLUTION TO A 

QUADRATIC EQUATION 

Let us demonstrate how to compute observable quantities within the field theory using the simplest example of the 
mean density of states. It is given by 



p(E) = (-2m)- 1 ((G« n (E)) - (Gt(E))) = ± ( V* {<p?tf + e~ m . 



(157) 



One can check that the pre-exponent in Eq. (|157p can be represented as 



(^V„ + = STr[n* n ® *„], 



(158) 



where 



n = u R - rr\ 



a o o o\ 





vO 0/ 



IT 



/0 ON 



10 

Vo o o 0/ 



(159) 
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Now we introduce an infinitesimal source field h n and add to the action Eq. (|140[) a term 

SF[Q,h] = -i^/iiSTrpQ,]. 

One can easily check that the density of states is given by a differentiation of the partition function with respect to 
the field h: 



p ( E ) = ^-.4- I™ e- F ^ 
Fy ' 2m dh n J 



, F[Q,h} = F[Q}+SF[Q,h}. (160) 

h— >0 



Note that the additional term in the action proportional to h enters exactly like the term proportional to w, so that 
in the final action of the sigma-model Eq. (1152p one can simply substitute 

|A->|a + /i„IL (161) 

Then Eq.(fT60| results in: 

p(E) = p (E) \JvQ STr[±IQ n ] e~ F ^\ (162) 

where po(E) is given by Eq. (|153[) and F[Q] is given by Eq. (|152|) at u> = 0. We see that the quantity po(E) which 
appear in Eq. (|152|) from the solution of the quadratic saddle-point equation Eq. (|148|) is not accidentally of the form 
of a semi-circle as the mean density of states is proportional to it. In the case of the Wigner-Dyson ensemble the 
functional F[Q] at ui = is simply zero and the integral in Eq. (|162[) is a constant independent of E. Thus we conclude 
that the semicircle law appears in this formalism from the solution of a quadratic saddle-point equation. 

Other statistics such as the two-point correlation functions can also be easily computed using the formalism of the 
nonlinear sigma model, however not so simply as the semi-circle law. For this one needs a proper parametrization of 
the matrix Q which resolves the constraints Eqs. (|150[) . (|151[) . 



XV. SYMMETRY OF SUPER-MATRICES Q AND Q. 

Let us return back to the derivation of the functional representation in terms of Q. It appears 12 that by a change 
of variables: 

ip R ' A = ±i^e ±iv/2+ia {l-\x R/ AXR/A) (163) 



p > = ±i^A 1/2 e xr/a, 

where A1/2 > 0, < (p < 2ir, < fl < ir and Xr/a, X*r/ A are * ne new anti-commuting variables, one can represent 
Q = ^ (g) ^ in the following form: 

-1 _ fu R 0\ f^RR Z RA \ fu^ 1 



II 



A 



The beauty of this form is that the commuting and anti-commuting variables are separated by factorization. Namely, 
the outer matrices U, U' 1 containing 2x2 matrices u R / A and W^, A 



{-Xr/aXr/a -Xr/a \ ( l ~ \Xr/ A Xr/a X* r / A \ 



ur/a = ( I , u R ) A = I (165) 

Xr/a 1 + 2X*r/aXr/a/ bf \ -Xr/a 1 + ^Xr/aXr/a) 



are made of the anti-commuting variables. The inner matrix E 



E 



^RR ^RA\ f^BB _0 



t AR t AA ) RA \ t FF/B ,. 
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which is diagonal in the FB space, contains only commuting variables with only BB sector non-zero: 



2 



RA 



One can show that the factorized form Eq. (|164p is common to both the field Q and the dual field Q, with matrices 
U, U" 1 being exactly the same. However the structure of the inner matrices E and E are different. Efetov has shownii 
that the constraints Eqs. (|150[) . (|151[) give rise to the following structure of E: 



J BB 



A VA 2 - 1 e iv 

-VA 2 - 1 e~*P -A 



(167) 



RA 



X F a/1 - A 



ltp F 



RA 



where A > 1, — 1 < Xp < 1, and ip, ipF G (0, 2n). 

To make practical calculations possible we also give (without derivation) the expressions for the Jacobians of the 
transformation from original variables to the variables of the above paramerization. They are 

J ^ = TYT"' ( 169 ) 
4 MM 

and 

J M = 8(A^F' (170) 

for the theories with coupling matrices A" 1 and Ay, respectively. 

The matrices Ebb of the structure Eq. (|166p as well as the matrices Ebb of the structure Eq. (|167[) can be diago- 
nalized by the pseudo-unitary rotation R: 

Ebb = RDR-\ E B s = RDR' 1 , R- 1 =A 2 R t A 2 , (171) 

where 

M |Al ~ A2l o Al ~ A %.-A 2 iV -*>)■ (i- .)- 

It is clear that the rotation matrix i? can be multiplied by a diagonal matrix 

e v 



^ € U(l) ® t/(l) 



without violating the condition of pseudo-unitarity and without changing the matrix E bb or Ebb ■ To eliminate the 
redundant degrees of freedom (which lead to the divergency of the functional integrals) the group of pseudo-unitary 
matrices {7(1, 1) should be factorized as R (U(l) ® U(l)), where R being a factor-group [/ff^'^i) ■ 
On top of that the diagonal matrix D has a free parameter 

Ai — A2 £ R.. 

The complete symmetry of the manifold of matrices Ebb and thus the complete symmetry of Q is: 

- U(l,l) 



U(1)®U(1) 

(7(1)®(7(1) ' 



In contrast to that the symmetry of matrices Ebb is simply Ebb € TfTT^kfjn\- Its counterpart E^p has the symmetry 



Tiff G 1^(^7(1) as ^ can ^ e di a g° nanze d by the unitary rotation matrix Rp. The complete symmetry of the Q field 
in the Efetov's nonlinear sigma-model is therefore: 

U(l,l) U(2) 
Qe U(1)®U(1) ^ U(1)®U(1)' 

Note by passing that the number of independent variables in Q and Q is different. While both have 4 anti-commuting 
variables, the number of commuting variables is 2 + 2 = 4 for the filed Q and 2 + 1 = 3 or the field Q. 

Thus we see that the duality transformation and the saddle-point approximation not only invert the coupling matrix 
Ay but also change the symmetry of the target space. Such type of duality is encountered in the string theory and is 
called T — duality. 
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XVI. EIGENFUNCTION STATISTICS 



In this section we show how to compute eigenfunction statistics using the field theory formalism. As usual, the 
starting point is to express the physical quantity of interest in terms of the Green's functions. To this end we study 
the product: 

^-[^[^-(e^b)' (E^)". <>«> 

where we used the representation of Green's functions in terms of exact eigenfunction ^(n) and exact eigenvalues 
E n of a random matrix Hamiltonian. Let us multiply Eq. (|173p by an infinitesimal S l + m ~ 1 average over realizations 
of the random matrix ensemble and do the limit 5 — > 0. This trick singles out only one state of the double sum, the 
one that is accidentally at the energy E: 

hm 5^ Kl m = lim /V |^(„)|^)^-i \ 

The smallness of the interval \E — Ei\ ~ 8 is the reason why the power of 8 is I + rn — 1 and not I + m. Indeed, let 
the joint probability distribution function for *i(n) and Ei be P(^i,Ei). It is a smooth function of E n which does 
not change at a scale 8^0. Then averaging in Eq. (|174|) can be performed as follows: 



E 



/I \T/ \2(l+m) r 



where 

r+ca 

(E - Ei + i8) l (E - Ei - iS) m ~ K "" J "" (J-l)l(m-l)! 

Let us define also the moment of the | S&i (rz) | 2 at an energy E: 



r [ + °° dEi , 9 «i-( I+m ) 9 . m _, q + m-2)! 



(I^| 2P )^ = ^) ^ W^-^)/ = E / / ^m,^)l^| 2p ^-^i)- (177) 

Comparing Eqs. (|177[) . (|175[) we arrive at: 

(l^i 2( ' +m) ) g = i l - m ~ 2) ] )l 5 hm {(2^) ^ +™- 1 ([Gl(i^; + ^<5)] ^ [Gl(^;-^<5)]™)}. (178) 

This is the expression of the eigenfunction moments in terms of the retarded and advanced Green's functions we were 
looking for. One can see that any non-trivial moment m + / > 1 requires a non-trivial limiting procedure. 

The next standard step is to represent the average of the Green's functions in terms of the functional integral. It 
begins with the standard representation similar to Eq. (|157[) : 

i-l r 

n R *,n R Y (m A *,n A \ m P ~ F \Q] 



[GZ n (E + l S)Y[G A n (E- l S)r = — J V#(<pg><pt<y(<p?rire-*W. (179) 

Then the analogy with Eq. (fT57|) would suggest that we write (p%*cp% = STr[H R Q], raise the hSTr[H R Q] into the 
exponent with the help of the ^-times differentiation with respect to the background field h and then switch to the 
super-matrix filed Q as in Eq. (|162j) . However, in trying to do these "standard" steps we make a mistake. The reason is 
that the field Q n is not slow-varying with n and the background field h n should also contain fast space variations. This 
is what makes a difference compared to the case of the constant in space symmetry breaking field ^wA in Eg. (|161[ > . 

One possible remedy^ is to single out the bi-linear combinations of cp* and (p n which do not contain fast space 
variations. We show how to do this for the product <p R *ip R tp A *(p A . As the result of averaging over random matrix 
ensemble should not depend on n (translational invariance on the average) one can do the sum over n and then divide 
the result by N. Switching to the Fourier-components ip(p)we can represent this sum as 

£ V R *{pi)v R {-px + q)v A *( P 2)<p A (-P2 - q) = E ^(Pl)^ B (-P2 + q)<p A *(P2)<p A (-Pl - ?)• 
Pi,P2,q pi,P2,q 
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Two sums in the above expression is a mere re-labeling of momenta, all what is really needed is that the sum of all 
momenta is zero. However, this re-labeling becomes a non-trivial operation if one assumes that the momentum q is 
small. In assuming so we select a definite domain of summation such that the corresponding bi-linear combination of 
if is slow varying in space. Then one single sum can be presented as 

^ tp R *( Pl )tp R ( P2 )<p A *(p 3 )cp A ( P4 ) S( Pl +P2+P3+PA) = (180) 



¥> R *(pi)<p R (-Pl + q)¥> A *(p2)v A (-p2-q)+ V R *{Pi)v R {-P2 + q)v A *{P2)V A {-Pi-<l)+ remainder. 

P1,P2,<?<1 PltJ>2>«<l 

In the first term of Eq. (|180p the bi-linear combinations if R * ( P i)ip R (— p\ + q) and tp A * ( P 2)<p A (— pi — q) are slow, while 
in the second term slow are the combinations tp R * (pi)(p A (— p\ — q) and <p A *(p2)<P R (—P2 + <?)• In the remainder we 
collect all terms where there is no bi-linear slow combinations. The meaning of the above procedure of singling out 
the slow bi-linear combinations is that only such combinations lead to the divergent functional integral in the limit 
when E + — E_ = 2i5 tends to zero. The average of the remainder is not singular and can be neglected. 

Performing this procedure in Eq. (|179[) one obtains (I + to)! possibilities to break the product (<fi R *<f R ) 1 (f A *f A ) m 
into the product of slow bi-linear combinations. All of them appear to make the same contribution to Eq. p79|) . Thus 
one can consider only one such term, do all the standard manipulations with the source fields as we explained above 
for the case of the mean density of states and multiply the result by q\. The final result for the simplest choice m = 1 
is: 

(\*n\ 2k )E = -~ lim|(27rp5) fc - 1 J VQ (STrp^QJ)^ 1 STr[U A Q n ] e~ F ^\ . (181) 

One can see that it is k\ times larger than the one obtained by the "naive" manipulations with the background field. 
We spent some time to go into detail of this subtlety in order to show that sometimes "exact" manipulations with the 
background fields are dangerous if the fast varying components of the fields are treated improperly or simply omitted. 

This is the result of a saddle-point approximation used in the derivation of the nonlinear sigma-model. No such 
danger appear for the dual representation which did not involve any approximation: 

(|*»| 2 *)b = - 2(fc 1 _ 1) , lim ^p5f- 1 J VQ (STrTT^])^ 1 STr[II A Q n ] e - f[<31 j . (182) 

One can do one more step without specifying the functionals F[Q] and F[Q] using the fact that the structure of 
dependence of the Q and Q fields on the anti-commuting variables Eqs. (|164p . (|165p is the same. To this end we 
definei^the generating functions Y[Q] as the functional integral of e - ^^ done over all the super-matrices Qi : except 
the one at a space point n: 

Y[Q n ] = I VQe- F ^. (183) 
If the generating function is known the eigenfunction moments are given by the integral over one single super-matrix 

Qn • 

One can show quite generally that this function does not depend on the anti-commuting variables. Then the 
integration of the anti-commuting variables is very simple as it involves only the pre-exponent in Eqs. (ll8ip . (|182p . As 
the result of this integration the additional factor [k — 1) appears in these equations. However, the main thing is to 
understand how it comes that the infinitesimal factor 5 k ~ 1 is compensated by the integral over the super-matrix Q n . 
There is only one scenario of for this to happen in the framework of the nonlinear sigma-model: this is to absorb 
S into the variable A — > 6X which can take arbitrary large values. For this the pre-exponent in Eq. (|18ip must be 
proportional to A fc ~ 2 in the limit of large A and also the generating function Y(Q n ) — Y(u) must be a function of one 
single variable u = 2npd\. One can show that this is indeed the case: 

(\* n \ 2k ) E = k -^A I™ duu k -*Y{u). (184) 







N 

This is a remarkable formula, as it implies that the distribution function of l^l^ for any unitary ensemble (3 = 2 is: 

(185) 



u=|*| ; 
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For the dual theory Eq. (|182[) . the generating function defined similar to Eq. (|183p may depend on the two 

variables. This is because there are not one but two non-compact variables Ai and A2 that may take arbitrary large 
values. Introducing new variables 

s = (Ai+A 2 ), r = (Ai-A 2 ), 

one obtains: 

1 1 poo r + 00 

(l^! 2 % = 4^7^)! J, d * L * * k - 2 ns,r). (186) 

Let us apply Eq. (|185p to the simplest case of the eigenfunction statistics in the Wigner-Dyson random matrix theory. 
In this case the variables of the super-matrix Qi are locked to their values at i = n. Thus there is no integration in 
Eq. (|183p whatsoever and we obtain: 



Y[Q] = exp [-irNpS STr[AQ n ]] -► exp(-27rp<5A N). (187) 
Then Eq. (|185p immediately gives for (3 = 2 Wigner-Dyson RME the Gaussian eigenfunction distribution: 

V(\V\ 2 ) = Ne- NW '. (188) 

Note that the Gaussian form of the distribution function is not a consequence of the Gaussian distribution of the 
entries of H but rather a consequence of the central limit theorem at any distribution of independently fluctuating 
entries which variance matrix Ay does not depend on i — j. One can show that for the orthogonal Gaussian ensemble 
P = 1 it deviates from the Gaussian: 



nm = ] j^ W e-^ (189) 

The simplest non-trivial application 12 of Eq. (|186|) for the problem that cannot be treated by the nonlinear sigma- 
model is calculating the eigenfunction distribution function for the one-dimensional Anderson model. This model is 
described by the random matrix Hamiltonian 

— Si, Hz.zil 1) 

where e$ is the Gaussian random variable with the variance w <C 1 (weak disorder case) . Outside the center of the 
band E = the result for the eigenfunction distribution function is amazingly simple: 

y. e -v loc |*| 2 

w 2 )4V' (190) 

where V\ oc is the localization radius. This distribution is not normalizable and should be cut a small values of l^l 2 . 
However, there is another way of normalizing it. This is the requirement that (|'I' I i| 2 ) = N . The first moment of 
the distribution is perfectly well defined and gives the above pre-factor. 
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